On the constraints violation in forward dynamics of multibody systems

It is known that the dynamic equations of motion for constrained mechanical multibody systems are frequently formulated using the Newton-Euler’s approach, which is augmented with the acceleration constraint equations. This formulation results in the establishment of a mixed set of partial differential and algebraic equations, which are solved in order to predict the dynamic behavior of general multibody systems. The classical resolution of the equations of motion is highly prone to constraints violation because the position and velocity constraint equations are not fulfilled. In this work, a general and comprehensive methodology to eliminate the constraints violation at the position and velocity levels is offered. The basic idea of the described approach is to add corrective terms to the position and velocity vectors with the intent to satisfy the corresponding kinematic constraint equations. These corrective terms are evaluated as function of the Moore-Penrose generalized inverse of the Jacobian matrix and of the kinematic constraint equations. The described methodology is embedded in the standard method to solve the equations of motion based on the technique of Lagrange multipliers. Finally, the effectiveness of the described methodology is demonstrated through the dynamic modeling and simulation of different planar and spatial multibody systems. The outcomes in terms of constraints violation at the position and velocity levels, conservation of the total energy and computational efficiency are analyzed and compared with those obtained with the standard Lagrange multipliers method, the Baumgarte stabilization method, the augmented Lagrangian formulation, the index-1 augmented Lagrangian and the coordinate partitioning method.The first author expresses his gratitude to the Portuguese Foundation for Science and Technology through the PhD grant (PD/BD/114154/2016). This work has been supported by the Portuguese Foundation for Science and Technology with the reference project UID/EEA/04436/2013, by FEDER funds through the COMPETE 2020 – Programa Operacional Competitividade e Internacionalização (POCI) with the reference project POCI-01-0145-FEDER-006941.info:eu-repo/semantics/publishedVersio

A study on wear evaluation of railway wheels based on multibody dynamics and wear computation

The wear evolution of railway wheels is a very important issue in railway engineering. In the past, the reprofiling intervals of railway vehicle steel wheels have been scheduled according to designers' experience. Today, more reliable and accurate tools in predicting wheel wear evolution and wheelset lifetime can be used in order to achieve economical and safety benefits. In this work, a computational tool that is able to predict the evolution of the wheel profiles for a given railway system, as a function of the distance run, is presented. The strategy adopted consists of using a commercial multibody software to study the railway dynamic problem and a purpose-built code for managing its pre- and post-processing data in order to compute the wear. The tool is applied here to realistic operation scenarios in order to assess the effect of some service conditions on the wheel wear progression

A parametric study on the dynamic response of planar multibody systems with multiple clearance joints

A general methodology for dynamic modeling and analysis of multibody systems with multiple clearance joints is presented and discussed in this paper. The joint components that constitute a real joint are modeled as colliding bodies, being their behavior influenced by geometric and physical properties of the contacting surfaces. A continuous contact force model, based on the elastic Hertz theory together with a dissipative term, is used to evaluate the intra-joint contact forces. Furthermore, the incorporation of the friction phenomenon, based on the classical Coulomb’s friction law, is also discussed. The suitable contact-impact force models are embedded into the dynamics of multibody systems methodologies. An elementary mechanical system is used to demonstrate the accuracy and efficiency of the presented approach, and to discuss the main assumptions and procedures adopted. Different test scenarios are considered with the purpose of performing a parametric study for quantifying the influence of the clearance size, input crank speed and number of clearance joints on the dynamic response of multibody systems with multiple clearance joints. Additionally, the total computation time consumed in each simulation is evaluated in order to test the computational accuracy and efficiency of the presented approach. From the main results obtained in this study, it can be drawn that clearance size and the operating conditions play a crucial role in predicting accurately the dynamic responses of multibody systems.Fundação para a Ciência e a Tecnologia (FCT

Spatial rigid-multi-body systems with lubricated spherical clearance joints : modeling and simulation

The dynamic modeling and simulation of spatial rigid-multi-body systems with lubricated spherical joints is the main purpose of the present work. This issue is of paramount importance in the analysis and design of realistic multibody mechanical systems undergoing spatial motion. When the spherical clearance joint is modeled as dry contact; i.e., when there is no lubricant between the mechanical elements which constitute the joint, a body-to-body (typically metal-to-metal) contact takes place. The joint reaction forces in this case are evaluated through a Hertzian-based contact law. A hysteretic damping factor is included in the dry contact force model to account for the energy dissipation during the contact process. The presence of a fluid lubricant avoids the direct metal-to-metal contact. In this situation, the squeeze film action, due to the relative approaching motion between the mechanical joint elements, is considered utilizing the lubrication theory associated with the spherical bearings. In both cases, the intra-joint reaction forces are evaluated as functions of the geometrical, kinematical and physical characteristics of the spherical joint. These forces are then incorporated into a standard formulation of the system’s governing equations of motion as generalized external forces. A spatial four bar mechanism that includes a spherical clearance joint is considered here as example. The computational simulations are carried out with and without the fluid lubricant, and the results are compared with those obtained when the system is modeled with perfect joints only. From the general results it is observed that the system’s performance with lubricant effect presents fewer peaks in the kinematic and dynamic outputs, when compared with those from the dry contact joint model.Fundação para a Ciência e a Tecnologia (FCT

Lubricated revolute joints in rigid multibody systems

The main purpose of this work is to present a general methodology for modeling lubricated revolute joints in constrained rigid multibody systems. In the dynamic analysis of journal-bearings, the hydrodynamic forces, which include both squeeze and wedge effects, generated by the lubricant fluid, oppose the journal motion. The hydrodynamic forces are obtained by integrating the pressure distribution evaluated with the aid of Reynolds’ equation, written for the dynamic regime. The hydrodynamic forces built up by the lubricant fluid are evaluated from the system state variables and included into the equations of motion of the multibody system. Numerical examples are presented in order to demonstrate the use of the methodologies and procedures described in this work.Fundação para a Ciência e a Tecnologia (FCT

Comparison and Implementation of a Rigid and a Flexible Multibody Planetary Gearbox Model

We propose algorithms for developing (1) a rigid (constrained) and (2) a flexible planetary gearbox model. The two methods are compared against each other and advantages/disadvantages of each method are discussed. The rigid model (1) has gear tooth reaction forces expressed by Lagrange multipliers. The flexible approach (2) is being compared with the gear tooth forces from the rigid approach, first without damping and second the influence of damping is examined. Variable stiffness as a function of base circle arc length is implemented in the flexible approach such that it handles the realistic switch between one and two gear teeth in mesh. The final results are from modelling the planetary gearbox in a 500 kW wind turbine which we also described in Jørgensen et.al (2013)

Modeling of the condyle elements within a biomechanical knee model

The development of a computational multibody knee model able to capture some of the fundamental properties of the human knee articulation is presented. This desideratum is reached by including the kinetics of the real knee articulation. The research question is whether an accurate modeling of the condyle contact in the knee will lead to reproduction of the complex combination of flexion/extension, abduction/adduction and tibial rotation ob-served in the real knee? The model is composed by two anatomic segments, the tibia and the femur, whose characteristics are functions of the geometric and anatomic properties of the real bones. The biomechanical model characterization is developed under the framework of multibody systems methodologies using Cartesian coordinates. The type of approach used in the proposed knee model is the joint surface contact conditions between ellipsoids, represent-ing the two femoral condyles, and points, representing the tibial plateau and the menisci. These elements are closely fitted to the actual knee geometry. This task is undertaken by con-sidering a parameter optimization process to replicate experimental data published in the lit-erature, namely that by Lafortune and his co-workers in 1992. Then, kinematic data in the form of flexion/extension patterns are imposed on the model corresponding to the stance phase of the human gait. From the results obtained, by performing several computational simulations, it can be observed that the knee model approximates the average secondary mo-tion patterns observed in the literature. Because the literature reports considerable inter-individual differences in the secondary motion patterns, the knee model presented here is also used to check whether it is possible to reproduce the observed differences with reasonable variations of bone shape parameters. This task is accomplished by a parameter study, in which the main variables that define the geometry of condyles are taken into account. It was observed that the data reveal a difference in secondary kinematics of the knee in flexion ver-sus extension. The likely explanation for this fact is the elastic component of the secondary motions created by the combination of joint forces and soft tissue deformations. The proposed knee model is, therefore, used to investigate whether this observed behavior can be explained by reasonable elastic deformations of the points representing the menisci in the model.Fundação para a Ciência e a Tecnologia (FCT) - PROPAFE – Design and Development of a Patello-Femoral Prosthesis (PTDC/EME-PME/67687/2006), DACHOR - Multibody Dynamics and Control of Hybrid Active Orthoses MIT-Pt/BSHHMS/0042/2008, BIOJOINTS - Development of advanced biological joint models for human locomotion biomechanics (PTDC/EME-PME/099764/2008)

Development of a planar multi-body model of the human knee joint

The aim of this work is to develop a dynamic model for the biological human knee joint. The model is formulated in the framework of multibody systems methodologies, as a system of two bodies, the femur and the tibia. For the purpose of describing the formulation, the relative motion of the tibia with respect to the femur is considered. Due to their higher stiffness compared to that of the articular cartilages, the femur and tibia are considered as rigid bodies. The femur and tibia cartilages are considered to be deformable structures with specific material characteristics. The rotation and gliding motions of the tibia relative to the femur can not be modeled with any conventional kinematic joint, but rather in terms of the action of the knee ligaments and potential contact between the bones. Based on medical imaging techniques, the femur and tibia profiles in the sagittal plane are extracted and used to define the interface geometric conditions for contact. When a contact is detected, a continuous non-linear contact force law is applied which calculates the contact forces developed at the interface as a function of the relative indentation between the two bodies. The four basic cruciate and collateral ligaments present in the knee are also taken into account in the proposed knee joint model, which are modeled as non-linear elastic springs. The forces produced in the ligaments, together with the contact forces, are introduced into the system’s equations of motion as external forces. In addition, an external force is applied on the center of mass of the tibia, in order to actuate the system mimicking a normal gait motion. Finally, numerical results obtained from computational simulations are used to address the assumptions and procedures adopted in this study.Fundação para a Ciência e a Tecnologia (FCT